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We provide an introduction to the use of ion crystals in a Penning trap P [2 [S] g] for 
experiments in quantum information. Macroscopic Penning traps allow for the contain- 
ment of a few to a few million atomic ions whose internal states may be used in quantum 
information experiments. Ions are laser Doppler cooled [l], and the mutual Coulomb re- 
pulsion of the ions leads to the formation of crystalline arrays [5][6l[3[8]. The structure 
and dimensionality of the resulting ion crystals may be tuned using a combination of 
control laser beams and external potentials [9] 110) . 

We discuss the use of two-dimensional ^Be"*" ion crystals for experimental tests of 
quantum control techniques. Our primary qubit is the 124 GHz ground-state electron 
spin flip transition, which we drive using microwaves 1111 1121 . An ion crystal represents 
a spatial ensemble of qubits, but the effects of inhomogeneities across a typical crystal 
are small, and as such we treat the ensemble as a single effective spin. We are able to 
initialize the qubits in a simple state and perform a projective measurement [l] on the 
system. 

We demonstrate full control of the qubit Bloch vector, performing arbitrary high- 
fidelity rotations (tjt ~200 lis). Randomized Benchmarking |13| demonstrates an error 
per gate (a Pauli-randomized 7r/2 and tt pulse pair) of 8 ± 1 X 10~*. Ramsey inter- 
ferometry and spin-locking |14| measurements are used to elucidate the limits of qubit 
coherence in the system, yielding a typical free-induction decay coherence time of T2 ~2 
ms, and a limiting Tip ~688 ms. These experimental specifications make ion crystals 
in a Penning trap ideal candidates for novel experiments in quantum control. As such, 
we briefly describe recent efforts aimed at studying the error-suppressing capabilities 
of dynamical decoupling pulse sequences [l5l [Tel [TtI [Tsl [Tol [20l [2TI [22. .23, ,24l [25l |26] . 
demonstrating an ability to extend qubit coherence and suppress phase errors |lll 1121 . 
We conclude with a discussion of future avenues for experimental exploration, including 
the use of additional nuclear-spin-flip transitions for effective multiqubit protocols, and 
the potential for Coulomb crystals to form a useful testbed for studies of large-scale 
entanglement. 
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1 Introduction 

Trapped ions have proven to be among the leading technologies for the development of quan- 
tum information systems. Typically, experiments are performed using Paul traps which pro- 
vide confinement to charged particles via a ponderomotive potential arising from an applied 
radiofrequency field. A technological alternative to the Paul trap is the Penning trap which 
provides confinement via the application of static electric and magnetic fields [27] . These traps 
have served as the fundamental components of experimental studies in fundamental particle 
physics ^28^, ^29^, '30', '3T] , strongly coupled plasmas (relevant to astrophysics) [32l [33l [2] , mass 
spectroscopy [34 , ,35 j . and the engineering of precision frequency standards |36| . Penning traps, 
however, have received relatively little attention as an experimental apparatus of choice in 
quantum information experiments j37] . 

In this article we describe our ability to implement quantum control techniques with 
high fidelity using arrays of ^Be+ ions confined in a macroscopic Penning trap [11 El |4j. 
We describe the fundamental operational principles and techniques for ion confinement and 
control in Penning traps '27], and describe their applicability for experiments in quantum 
information. Our experiments employ two-dimensional planar crystals [6] of ~1000 ^Be+ ions 
in a magnetic field of ~4.5 T. At this large magnetic field we use a ground-state electron-spin- 
flip transition with a splitting of ~124 GHz as a qubit basis, and implement quantum control 
via a microwave system [Tl] [12] . We demonstrate through measurements of qubit coherence 
and characterization of the experimental system that the effects of inhomogeneities are small, 
and as such we treat the spin ensemble as a single effective qubit. 

We present a detailed characterization of measurement fidelity, qubit coherence, and the 
operational performance of Pauli and Clifford rotations in this system. Our measurements 
are based on demonstrations [38] of high-contrast Rabi flopping, controlled Larmor preces- 
sion, randomized Clifford benchmarking [13], Ramsey free-induction decay, spin-locking [l4] . 
and multipulse dynamical decoupling pSl|16l[iail8lll9l|20l|2ll|2a|23l21l25l|26]. Experi- 
ments reveal long coherence times and high operational fidelity, despite the sensitivity of our 
qubit transition to magnetic field fluctuations. In addition to the dominant electron spin-flip 
transition, we also demonstrate the capability of shelving quantum information in nuclear 
spin states with ~300 MHz level splitting. Our control system permits the demonstration of 
environmental bath engineering [Til IE] , enabling our qubits to mimic the dynamics of qubit 
technologies which experience markedly different classical noise environments. In summary, 
our system may be thought of as a model quantum system, or a precise spin resonance system, 
providing the ability to perform high-fldelity state initialization and projective measurement. 

The remainder of this article is organized as follows. We begin with an introduction to the 
basic operating principles of ion storage in a Penning trap in section |2l We then move on to 
a detailed description of our experimental apparatus, including our Penning trap, microwave 
system, and qubit manifold of choice in section |3l Descriptions of high-fldelity quantum 
control, and detailed system characterizations are provided in sections |4l|5] and HI followed by 
a brief discussion of future experiments in section [71 

2 Penning Trap Basics 

In this section we provide a brief overview of the functional basis for the Penning trap, and 
highlight some of its salient operational characteristics, largely following the presentation in 
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[27j . The Penning trap uses static electric and magnetic fields in order to provide confinement 
for charged particles; an axial, constant magnetic field oriented along the Z-axis, B — BqZ, 
and an electric field that provides a harmonic potential well at a locus along this axis. The ideal 
trap consists of a central ring electrode that is a hyperbola of revolution and two hyperbolic 
endcaps, with distance from the center of the trap to the ring tq, and distance to the endcap 
zq. The electrode surfaces are constructed to be equipotential surfaces of the electric potential, 

$ ^ (2z^ -x^- y^) , (1) 

2(z2+r2/2)^ ^' 

where Uq is the fixed voltage between the ring and end-cap electrodes. Alternatively, simple 
cylindrical electrodes may be employed such that the potential is well approximated by Eq. [T] 
for distances small relative to the electrode dimensions. 

A charged particle is confined axially by the applied electric fields (Fig. [Iji) and radially 
by the magnetic field. Writing the particle's equations of motion in three dimensions wc find 
that the particle is axially confined in a harmonic potential with frequency 



UJ 



2 ^ 2qUo 



(2) 



The radial motion is described by a superposition of two distinct circular motions with fre- 
quencies, 

UJ+ ^ - {uJc + UJi) , W_ = - (cJc - Wl) . (3) 

Here, ujc — qBo/m is the cyclotron frequency for charge q and mass to, uj^ is known as the 
modified cyclotron frequency, cj- is the magnetron frequency, and lui = uj^ — 2a;^. Stable 
confinement is achieved when 

> (4) 

Thus we see that a charged particle in a Penning trap is harmonically confined in the 
axial direction, but in the radial plane the particle undergoes complex motion, and traces 
out what is known as an epitrochoid, due to the high-frequency cyclotron motion and the 
slower magnetron motion. The x — y motion can be periodic or quasiperiodic depending upon 
whether or not the ratio cj_|_ /u- is rational. 

3 The NIST Penning Trap 

The NIST Penning trap consists of a stack of cylindrical electrodes as pictured in Fig. [Hand 
described in detail in [l] [3] . Charged particles are confined at the symmetry point of the 
stack along the axial direction, as shown in the inset to Fig. [T^. Due to E x B drift from the 
axial magnetic and radial electric fields, the cloud of charged particles rotates about its axes 
at a frequency greater than tj_ . Segmented electrodes (Fig. [1^) are used to apply a rotating 
dipole or quadrupole potential know as a "rotating wall." Using this external potential it is 
possible to phase-lock the rotation of trapped ions to the wall, adding stability to the system. 

The entire trap structure and enshrouding vacuum envelope is inserted into a 4.5 T room- 
temperature-bore superconducting magnet with a vertical field orientation. Laser access par- 
allel and perpendicular to the magnetic field axis is provided using apertures and mirrors 
mounted on the exterior of the glass vacuum envelope. Similarly, optical access for imaging 
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is provided both along the field axis ( "Top- view" ) and perpendicular to it ( "Side-view" ) . For 
the quantum information experiments described in this article, a side-view imaging system 
using either a spatially resolving camera or a phototube were employed to collect scattered 
fluorescence from ^Be"*" ions as state-selective readout, as will be described in the next section. 
Top-view imaging was also employed in the characterization of ion crystals as will be reviewed 
in the next section. 

3.1 Trapping and Controlling ^Be^ Ion Crystals 

Previous experimental studies have demonstrated that it is possible to ionize, confine, and 
laser-cool up to ~10^ ^Be+ atoms in our Penning trap. Originally, these experiments focused 
on microwave frequency standards and the physics of cold, one-component plasmas [39l [1] , 
appropriate in certain circumstances for studies of dense astrophysical objects. We present 
here an overview of the salient physics of these systems that will assist in our understanding 
of the Penning trap as an experimental testbed for quantum control. 

When Doppler laser cooled to sufficiently low temperatures, trapped ^Be+ ions crystallize 
due to their mutual Coulomb interaction [Sj [6l [3] . This phenomenon is predicted in a model 
for a classical one-component plasma, and is governed by the Coulomb coupling constant [40j , 



47reo (3/47m)^ /cbT' 

where is the permittivity of free space, e is an ion's charge, fc^ is Boltzmann's constant, 
T is the ion temperature (which does not include the effects of ion rotation in the confining 
potential), and n is the plasma density. The plasma density is set approximately by the phys- 
ical extent of the ion cloud as limited by the confining potential. Reducing the temperature 
relative to the strength of the Coulomb interaction increases P, and above P w 170 theory 
predicts that Coulomb repulsion dominates thermal motion and leads to the formation of a 
body-centered-cubic (BCC) three-dimensional array [3]. 

The NIST Penning trap yields an ion cyclotron frequency o^c ~ Stt x 7.61 MHz at 4.46 
T, an axial frequency uj^ ~ 2Tr x 799 kHz and a magnetron frequency uj^ = 2t: x 42.2 kHz 
at a trap voltage of 1 kV. We have previously demonstrated the transition to BCC order in 
Doppler cooled ^Be+ plasmas by both direct imaging [6l |3] and Bragg scattering 5 studies 
using scattered light from the cooling laser. Calculations indicate that we achieve a coupling 
constant P « 2000, for an ion density n — A x 10^ cm~'^, and ion temperature T « 1 mK. 
While similar crystallization has been observed in RF Paul traps [41] , limitations imposed by 
the onset of RF heating from trapping fields have prevented the realization of large crystalline 
arrays with stable BCC ordering. 

In thermal equilibrium the ion plasma rotates rigidly in the confining potential. The 
rotational frequency is controllable via both an external rotating potential iQiilOj and a cooling 
laser beam applied perpendicular [T] to B. For example, the perpendicular laser induces a 
torque on the plasma due to radiation pressure, which can be used to increase or decrease the 
rotational frequency. The rotation frequency, w^, controls the shape and dimensionality of the 
ion plasma. This dependence arises from the fact that the radial confinement associated with 
ion motion in the magnetic field depends on the rotational velocity [TJ U] . Thus by changing 
LUr it is possible to adjust the effective radial trapping potential. The ion plasma adjusts 



(5) 
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Fig. 1. NIST Penning trap, (a) Photograph of trap electrode structure showing loading trap 
zone (bottom) and experimental zone (top). Segmented electrodes in the top half of the image 
are used for application of a rotating-wall potential. (Inset) Schematic of experimental zone 
electrode layout, confining electric fields, ion location (cloud at center), and principal cooling laser 
direction. (b)Schematic of system providing details of laser access/orientation, imaging system, 
and microwave system. Top-view imaging system not shown, (c-d) Top and side view images 
of resonant fluorescence from crystallized ^Be"*" ion plasmas. Top view shows hexagonal order, 
acquired stroboscopically, by gating on the rotating-wall potential. Multiple ion planes are present, 
leading to differing degrees of focus among imaged ions. Side view shows a different crystal with 
a single plane of ions. 
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its physical extent to the strength of the radial confinement: for weak confinement near the 
lower bound cUr w uj^ the cloud assumes an oblate form, while for uj « lUc/'2 the cloud spins 
up to an oblong cigar-like shape along the magnetic field axis. As the rotational frequency 
increases up to tOr ~ w+, the confining centripetal acceleration becomes insufficiently strong 
to maintain the axially extended shape of the ions and the cloud relaxes towards an oblate 
shape. At higher rotational frequencies the balance of the outwardly directed radial electric 
force and the confining Lorentz force does not provide confinement and the radial extent of 
the cloud diverges. For the quantum information experiments described in this work we focus 
on oblate clouds at low rotation frequencies with ~f 00-fOOO ions in one or two planes, easily 
resolvable via side-view imaging (Fig.[TJi), and having radial extent of a few hundred microns. 

Once the rotational frequency is coarsely set via the perpendicular laser, the ions are 
locked to an external dipole "rotating wall" potential [9l [10] , providing long-term stability 
over many hours. Figure [1]; shows a top view image of an ion array obtained with a long 
integration time (~30 s). This image was obtained with an intensified CCD camera, strobing 
the camera synchronously with the applied rotating potential used to control the plasma 
rotation frequency. The figure is therefore an image of the ions in the frame of the rotating 
wall potential. Ions are well resolved in the array center, but not near the edges of the array. 
As discussed in |42| this difference is due to small slips of the planar crystal relative to the 
rotating wall potential. Studies of the "stick-slip" dynamics of these crystals |6l |42] have 
revealed that the ions experience abrupt reorientations of the crystal structure in the frame of 
the rotating wall, separated by stable phase-locked rotation for periods many tens of seconds 
long. 

The dominant ion loss mechanisms limiting the ion cloud lifetime are background collisions 
in which ^Be+ is converted to BeH+ or BeOH+ [4 . These molecular ions remain part of 
the crystallized plasma, and are sympathetically cooled by the ^Be+ ions. These so-called 
"heavy- mass" ions centrifugally separate from the lighter ^Be"^, and due to their modified 
level structures do not significantly participate in any quantum control experiments. Ion 
collisions with background gas molecules typically reduce the number of ^Be+ ions by half 
over a timescale of one to two days. 

3.2 ^Be+ Level Structure at 4.5 T 

The level structure of ^Bc+ ions is modified by the presence of an external confining magnetic 
field with strength ~4.5 T. The relevant energy levels are depicted schematically in Fig. [2] 
[1]. Our primary qubit states are defined by two states with 124 GHz separation within the 
2s^5'i/2 manifold, 

\F = 2, mj = +3/2, mj = +1/2) = |t) ^ |F = 2, m/ = +3/2, m., = -1/2) = \i) , (6) 

where m/ and mj are the nuclear and electron spin projections along the quantizing field, 
and F is the hyperfine quantum number. Around 4.5 T, the qubit states diverge linearly at 
28 MHz/mT, making this transition sensitive to magnetic field fluctuations. The "bright" |t) 
state is coupled to the 2p^P3/2(™j = +3/2) manifold by a cycling transition at 313 nm. This 
transition is employed for Doppler cooling (when red-detuned) , plasma control (via radiation 
pressure in a perpendicular conflguration) and state-selective readout, since |t) is optically 
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bright. The "dark" state ||) is minimally affected by the radiation on the cycling transition ^, 
but can be rapidly repumped by another transition near 3f 3 nm to the 2p2p3/2(TOj = +1/2) 
manifold. 



^Be^ @ B=4.5T 



P_2p F=o, 1,2,3...:::: 



2S Si 



/2 



F=2 





CD 


Q-! 


O 

C D3 


E; 


CD C 


3 ! 


o = 


o ' 


(n o 


Rei 


Fluore 
Co 



Qubit 



^+3/2.-1/2^ 



i) 



Fig. 2. Level structure of ^Be+ at 4.5 T. Highlighted levels correspond to electron-spin-flip tran- 
sition used as primary qubit basis. Numbers in parentheses represent (m/,mj), the nuclear and 
electron spin projections along the quantization axis. The upper qubit state is purely defined 
by these quantum numbers while other states in the hyperfine manifold contain a small amount 
of admixing of other (m/,mj) states {mp is a good quantum number). Repump and Fluores- 
cence/Cooling transitions are near 313 nm. Qubit splitting is near 124 GHz. 



3.3 Microwave and RF Control Systems 

Transitions between the qubit states are driven by a home-built microwave system [111 I12j as 
depicted in Fig. [3l The system is based on a Gunn diode oscillator which outputs approxi- 
mately 30 mW at 124 GHz. This oscillator is phase-locked to the sum of the eighth harmonic 
of a 15.5 GHz DRO and an external reference at ^77 MHz (for definitions of hardware ab- 
breviations see Fig.[3|). 

The DRO employs a 100 MHz quartz crystal oscillator as an ultra-stable, low-phase-noise 
reference. The bandwidths of the DRO PLL (to the 100 MHz crystal oscillator) and the 
124 GHz Gunn diode oscillator PLL (to the DRO) are 100 kHz and several times 100 kHz, 
respectively. Accordingly, the phase noise of the microwave source for frequency offsets less 
than 100 kHz is determined by the phase noise of the 100 MHz crystal oscillator. Our 100 MHz 
crystal oscillator has a maximum phase noise specification of -125 dBc/Hz, -150 dBc/Hz, -175 



^The dark state is repumped to the bright state by the ~313 nm cooling light on a timescale of 5-10 ms. 
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dBc/Hz, and -175 dBc/Hz at 100 Hz, 1 kHz, 10 kHz, and 100 kHz respectively. This implies 
a maximum integrated phase noise of ~107 dBc over the experimentally relevant interval 
[100 Hz, 100 kHz]. The phase noise of the microwaves is increased by 201og(A^) = 62 dB 
where A^=1240 is the step- up factor from the 100 MHz oscillator to the 124 GHz microwaves, 
leading to an integrated phase noise of ~ —45 dBc at the 124 GHz carrier over the same 
interval. 

The ~77 MHz reference is provided by a computer-controlled direct digital synthesizer 
(DDS) which is referenced to a 1 GHz clock derived from the 5 MHz signal of a passive 
hydrogen maser. Changes in the ~77 MHz signal are tracked linearly by the Gunn Diode, 
making this reference the primary means of tuning the output of the microwave system. The 
DDS output may be switched between four pre-programmed output profiles via TTL inputs, 
each with a user-defined frequency and phase, and the Gunn diode PLL tracks changes in 
DDS output within ^5 /is. 

The configuration described above is engineered for maximum amplitude, phase, and fre- 
quency stability. For quantum control experiments, however, it is occasionally desirable to 
engineer the noise environment to exhibit precisely defined spectral characteristics, as will 
be seen in section [5.41 As it is difficult to inject magnetic field noise into a superconducting 
magnet system, we instead engineer noise in the microwave control system |11[ 112] , which in 
the appropriate frame of reference is indistinguishable from the action of classical magnetic 
field fluctuations on the qubit states. We replace the ~77 MHz reference for the Gunn diode 
PLL with a frequency- modulated (FM) signal generated by a synthesizer (Fig.[3|D). A desired 
spectrum is generated numerically and Fourier transformed to produce a time-domain data 
set whose two-time correlation function reproduces the spectrum of interest. These data are 
sent to an arbitrary waveform generator whose voltage output is sent to the FM-input of 
the synthesizer, producing a carrier at 275 MHz with the desired noise characteristics. This 
carrier is mixed with a tunable ~197 MHz output from the computer-controlled DDS, band- 
pass filtered, and sent to the Gunn diode PLL, which linearly tracks the frequency-modulated 
reference to a bandwidth of ^100 kHz. 

Microwave power is output from the Gunn diode to a horn antenna via rigid waveguides. 
Power reaching the antenna may be tuned via a manual attenuator, and pulsed-control is 
afforded by a p-i-n diode switch that provides ~25 dB of isolation in the off-state. Using 
this hardware configuration, the only available pulse shape is approximately square, and 
the relatively long lengths of the relevant pulses (~ 100-200 /is) yield a nearly ideal square 
envelope. Output from the microwave horn propagates in free-space and is focused on the ion 
crystal using a custom teflon lens installed in the bore of the superconducting magnet. The 
microwave beam waist at the position of the ions is approximately 7 mm, allowing for large 
incident microwave power on the ions. The beam waist and microwave wavelength are large 
relative to the extent of the ion crystal, allowing for high driving-field homogeneity. 

An RF coil antenna is located outside of the vacuum envelope, in close proximity to the 
location of the ions, and is used to drive nuclear-spin-flip transitions as will be discussed in 
section 1^1 Coarse and fine synthesizer outputs are mixed and bandpass filtered to provide 
the appropriate carrier frequencies with resolution of ~1 mHz (Fig. [3J:). Both synthesizers 
are phase-locked to the output of a passive hydrogen maser, as used for the DDS. The carrier 
signal is amplified using a standard RF amplifier and power reaching the antenna is controlled 
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via a commercial RF switch under TTL control. 

Pulsed-control over the qubit drive and RF systems is achieved via a programmable-logic 
device. Details of the experimental control and pulse-sequencing systems, as well as the user 
interface, are provided in reference [l2j . 



(a) Microwave System (b) Noise Generation (cl RF System 




DDS/SYNTH 

Fig. 3. Schematic block diagrams of a) 124 GHz microwave system, b) Noise generation system, 
and c) RF system. PLL = Phase Lock Loop, DRO = Dielectric Resonator Oscillator, CPL = 
Directional Coupler, DDS = Direct Digital Synthesizer, ARB = Arbitrary waveform generator, 
SYNTH = frequency synthesizer, TTL = Transistor- Transistor Logic pulse. 



4 Quantum Control Experiments 

All experimental sequences begin with optical pumping ^ of the ions to the bright state via the 
cooling/detection and repump lasers near 313 nm. Pulsed exposure for ~100 ms, synchronous 
with line frequency initializes the ions to ||) with infidelity ^ 1 x 10~^, estimated from 
the scatter rate of the repumping laser and the |t) to ||) transition rate due to microwave 
leakage. Ion fluorescence is integrated over 50 ms using an //5 optical system on a phototube, 
an experimental microwave or RF pulse sequence is executed, and ion fluorescence is again 
measured for 50 ms. In all experiments, the ions are rotated to the dark state, at the end 
of the control sequence. The second measurement cycle is broken into five 10 ms bins and 
a linear fit to the count rate provides the fluorescence level at the end of the experimental 
control sequence. This technique accounts for repumping due to the detection laser from 
II) to It). The extracted fluorescence count rate is normalized to the bright-state count rate 
measured at the beginning of the experimental cycle to mitigate the effects of slow fluctuations 
in laser amplitude between experimental runs. Experimental data are generally presented as 
a normalized fluorescence count rate, which is commensurate with population in If), and is 
sometimes used as a measure of the qubit coherence, as will be described in the context of 
each experiment. Each data point presented in an experiment typically consists of an average 
of 20-50 individual experiments, all performed under approximately the same conditions. 
In the following subsections we will describe a series of quantum control experiments, 

^Optical pumping of the nuclear spins occurs due to mixing of different \mj, mj) states in the excited P-state 
manifold over a time ~ 1 s. For quantum control experiments using the 124 GHz electron-spin-flip transition 
this repumping does not need to be repeated because the nuclear spin does not leave the mj = -f 3/2 manifold. 
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primarily using the electron-spin-flip qubit basis. Each experiment consists of an appropriate 
microwave or RF pulse chain applied in the experimental sequence described above. These 
experiments will elucidate our quantum control capabilities, and will highlight the strengths 
and limitations of this experimental system. 

4-1 Driven Rotations 

We drive transitions between the electron-spin-flip qubit states using resonant microwave 
radiation focused on the ions as described in the previous section. In these experiments we 
apply a single microwave pulse with variable duration, tp, and measure the population in ||) 
as a function of pulse length. In Fig. [4^ we demonstrate Rabi flopping with contrast ^99.85 % 
(measured for tp = r^r) [431 144j with a 7r-pulse time, ~ 185 fis. The value of is tunable 
from roughly 90 /iS to 700 jis using a manually adjusted attenuator, and is extracted from 
fits to on-resonance Rabi-flopping curves. The measured contrast arises from a convolution 
of measurement and operational infidelities. Typical experimental conditions set « 200 /iS 
in order to simultaneously maximize flopping contrast and PLL stability. For shorter values 
of Ttt, reflections from the microwave switch mandate manual adjustment of a phase shifter 
(to mitigate reflections) on an hourly timescale. The phase of the DDS PLL reference may 
be changed phase-coherently by 90 degrees in order to affect ax and cry rotations, where 
(7; is a Pauli matrix. However, the absolute phase of the DDS reference is insignificant in a 
single-pulse experiment. 

The contrast of Rabi flopping curves can be seen to decay by ~l/e on a timescale of 
30-40 ms, as in Fig.|3)D. The form of the decay is not well-approximated by either a simple ex- 
ponential or a Gaussian, suggesting the loss in contrast is due to a complex interplay between 
multiple processes including long-term phase and amplitude instability of the microwave sys- 
tem ^44j. The integrated phase noise of the microwave carrier described earlier produces a 
phase deviation of a few tenths of a degree over a timescale of a few milliseconds, suggesting 
phase instability is significant, but not necessarily the only source of the observed decay in 
Rabi oscillations. As we will describe in Sec. [5l we believe that the frequency instability of 
the qubit transition due to fluctuations in the ambient magnetic field dominates expected 
instabilities in the microwave system. 

Rabi lineshapes in the frequency-domain are illustrated for 6 = w, 27r, 37r, and 47r rotations 
in Fig.[4j;. These data are obtained by applying a fixed-duration microwave pulse to the qubits, 
for various values of the ~77 MHz reference output by the DDS. Fits to these data incorporate 
the measured r^r and rotation angle, and a single free parameter, luq, the frequency of the 
DDS reference on resonance. The Rabi lineshape 22] is described by 



where lu is the DDS reference frequency. 

The resonance frequency drifts slowly due to flux leakage from our superconducting magnet 
at a relative rate of < 1 x 10~^ per day. Additionally, luq can occasionally be seen to shift 
suddenly by up to ^1 kHz, presumably due to the influence of external laboratory equipment 
in the building. Careful monitoring of ujq allows us to maintain long-term microwave frequency 
stability relative to the qubit transition of ^Ix 10^^ over any experimental measurement run. 
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Rabi Time (ms) nWave Ref. Frequency (Hz) 

Fig. 4. Driven qubit rotations using 124 GHz microwave system, a) Rabi flopping with = 185 fis. 
b) Decay of Rabi flopping contrast for long driving times. Fringe decay is not well approximated 
by either a simple exponential (dashed) or a Gaussian (dotted) decay, suggesting the possibility of 
a complex interplay of processes leading to loss of contrast . Inset) Graphical depiction of direction 
of driven rotations on the Bloch sphere, c) Rabi lineshapes for various rotation angles, determined 
by the rotation angle 9 (8 = n for ip = r^), as a function of the DDS reference frequency for the 
124 GHz microwave system. Solid lines are a fit to the Rabi lineshape (see text) with a single free 
parameter - the center frequency, uJo- 
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Ramsey Precession Time (ms) Ramsey Precession Time (ms) 

Fig. 5. Controlled Larmor precession, a) Pulse schematic for Ramsey experiments used to ob- 
serve Larmor precession. Microwaves (switch in off state) are detuned from resonance during the 
interpulse period by A, leading to precession at oj = 27rA. b-f) Qubit state population at the 
end of the Ramsey experiment as a function of precession time. Larmor precession frequency 
increases linearly with A as expected. Fringes show near unity contrast over measurement range. 
The observed 20 ^s offset is due to a delay period during which A = during free-precession, 
incorporated for stability of the PLL. Error bars shown in panel f) depict standard deviation of 
20 experiments averaged to produce each point. Standard deviation is minimized at the extrema 
of the Ramsey fringes, consistent with quantum projection noise. Note that actual error per data 
point is l/\/20 of the shown value. 
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4-2 Controlled Larmor Precession 

Rotations about the quantization axis in the qubit basis are achieved via detuning of the 
microwave drive frequency. We demonstrate controUed Larmor precession using a standard 
Ramsey interference experiment |44j . We begin a Ramsey experiment by rotating the Bloch 
vector of the qubits to the equatorial plane via a 7r/2x pulse. Here the subscript denotes the 
axis of rotation in the frame rotating at the unperturbed qubit splitting, = 124 GHz + loq. 
The microwave frequency is then detuned from resonance by an amount A, and the qubit 
is allowed to freely precess. At the end of the free-precession period another ■n/2x pulse is 
applied. The qubit state after the second it /2x pulse depends on the phase accumulated by 
the qubit in the rotating frame, which is manifested as a sinusoidal oscillation in the qubit 
state population in W) as a function of free-precession time. The period of the oscillation is 
A~^, and varying A modifies the oscillation period, as shown in Fig.O Setting the precession 
time to ^A~^ while detuned by A performs the Pauli gate nz- The Ramsey fringes measured 
for A = 100 kHz are plotted with error bars corresponding to the standard deviation of the N 
measurements averaged together for each data point (error bars are not normalized by ^/N). 
We observe that the standard deviation oscillates as a sinusoid with twice the frequency of 
the Ramsey fringes, and a 7r/2 phase shift. This yields maximum values of the standard 
deviation for states which have a 50 % probability of being projected to |t) or ||), and 
minimum deviations for states projecting with unit probability to either of the basis states. 
The observed periodicity and normalized values of the standard deviation are consistent with 
previous studies of quantum projection noise [45j . 

4-3 Randomized Gate Benchmarking 

We characterize the operational fidelity in our system using a simplification of the randomized 
benchmarking procedure introduced and demonstrated in reference |13j and subsequently 
employed in [46|, I47j . This procedure provides a simple mechanism to characterize the average 
fidelity of computationally relevant Clifford operators without the onus of performing full 
process tomography. Clifford operations implemented as driven 7r/2 gates about random axes 
are performed sequentially, with interspersed Pauli tt rotations performed about arbitrary 
axes 

in order to randomly shift the reference frame. A single computational gate in this con- 
struction is therefore composed of both a Clifford 7r/2 and a Pauli tt operation, each about 
a randomized axis, as depicted in the inset to Fig. [G] k runs of length I random computa- 
tional gates are performed, and the resulting fidelities are averaged for each value of I. For 
independent, depolarizing errors, the error probability for an Z-gate sequence, p;, increases 
exponentially with I. Fitting an exponential decay as a function of I to the average of the k 
randomizations thus provides a measure of the error per computational gate. 

We assemble a sequence of randomized computational gates and apply them to an initial 
qubit state |t). X operations are driven on resonance with r^r = 232.5 /is; driven Y operations 
are realized by shifting the phase of the ^77 MHz reference to the PLL provided by a DDS by 
90 degrees, using TTL inputs. Each operation is separated by a delay of 5 /us to ensure that 
the PLL tracks the phase changes of the DDS reference. Software keeps track of the expected 
outcome at the end of the computational sequence and applies an appropriate rotation to the 
Bloch vector such that in the absence of any error the state is rotated to \[) for measurement. 
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Computational Gate Number, / 

Fig. 6. Randomized benchmarking experiments using driven X and Y operations. Randomized 
sequence of I computational gates constructed as depicted in inset. Final rotation is deterministic 
and rotates qubit to dark state as described in the text. Small markers represent average fidelity of 
an individual randomized sequence of I gates, with each point consisting of 20 averaged measure- 
ments of the same sequence. Average fidelity defined as (1 — P\-\))- Large markers show average 
across k = 20 randomizations for each value of I. Solid line indicates weighted exponential fit to 
the data, and dashed lines indicate 68 % confidence level of the fit. (Inset) schematic depiction 
of the pulse sequence employed. All gates are separated by a 5 /iS delay during which microwave 
drive is turned off to allow the PLL to track phase jumps in the external reference signal 
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This is accomplished by generating an appropriate rotation matrix for each Chfford and 
Pauh operation, performing matrix multipUcation to simulate the string of computational 
gates, and generating an appropriate final rotation based on the expected state at the end 
of the computation. This implementation is an approximation to the prescription published 
in Ref. [13] which omits az operations as they simply require an arbitrary change of basis 
and can have much lower error rates than actual driven operations. Additionally, we perform 
all rotations with uniform sign. Despite these limitations this construction provides a good 
approximation to the procedure outlined previously, and focuses on the limiting operational 
fidelities associated with driven microwave transitions. 

In Fig. [6] we present results for fc = 20 random gate sequences with up to Z = 200 gates. 
We plot the average fidelity as a function of /; the fidelity decays to a value of 0.5 for complete 
depolarization. Fitting a simple exponential decay to the average fidelity as a function of I 
provides the solid line presented in the figure, yielding an average error per computational gate 
of 8± 1 X 10~*. Dashed lines indicate the 68 % confidence interval for the fit. These results are 
among the best reported to date for atomic systems, and the low error per computational step 
is largely due to the extraordinary stability of driven microwave rotations and the absence 
of spontaneous emission probability in microwave transitions, in contrast to laser-mediated 
gates. We believe the primary limitation on operational fidelity is due both to ambient 
magnetic field fluctuations and to small inhomogeneities in the magnetic field over the extent 
of the ion crystal, as will be discussed in Sec. 15.21 Additionally, pulse timing resolution of 50 
ns corresponds to a systematic 7r-pulse timing infidelity of up to '--^ 2 x 10^^. Driven 7r-pulse 
rotations, however, are sensitive to timing infidelity only to second order. 

These measurements are performed under typical experimental conditions, and no extra 
care has been taken to mitigate suspected sources of error. Previous experiments in supercon- 
ducting qubits [47] have demonstrated a dependence of pi on r^r, and, based on the relevant 
sources of microwave instability in our system, it should be possible to improve our extracted 
pi by reducing the length of driven rotations. However, in our current system, substantially 
reducing r^r has been observed to produce frequency instabilities in the microwave system at 
pulse edges. 

5 Qubit Coherence 

In this section we present a series of measurements designed to provide information about 
the limits of qubit coherence in our experimental system. We study free induction decay in 
a Ramsey interference experiment, relaxation in the rotating frame via spin-locking, and the 
suppression of dephasing by application of dynamical decoupling pulse sequences. In addition, 
we present an analysis of the effects of magnetic field inhomogeneities in our experimental 
apparatus. 

5.1 Free Induction Decay 

In addition to demonstrating controlled az operations, Ramsey free- induction decay mea- 
surements provide information about the coherence time of the qubit system [441 138j . In 
Fig. [7^ we show Ramsey fringes obtained using A — 100 kHz to 4 ms precession time. Fringe 
contrast decays as a function of precession time due to a loss of phase coherence between the 
microwave drive and the qubit system. 
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We consider the phase coherence of the qubit based on the following analysis. Starting 
with a qubit Bloch vector oriented along the Y-axis in the equatorial plane, we employ a 
measure of the coherence as W{t) — {ayit)) , where the brackets indicate an expectation 
value and the overline indicates an ensemble average [221 [211 [TTl [T2] . In a frame rotating at 
the qubit precession frequency, and in the absence of dephasing, the state remains oriented 
along the Y axis and W{t) ^ 1. In our system, classical magnetic field fluctuations form the 
dominant source of noise and hence provide the limiting mechanism for qubit coherence. 
The longitudinal relaxation time, Ti is effectively infinite in this system and will not be 
considered in our discussion. A Hamiltonian accounting for the presence of this noise may be 
written as [221 EH [111 [12] 

H=^[n + m]^z. (8) 

where VL is the unperturbed qubit splitting, and [3{t) is a time-dependent classical random 
variable describing the noise in the system. Under the influence of this Hamiltonian a state 
originally oriented along Y in the rotating frame may be written as 

\^{t)) = ^-'^ !>^'">^'' |t) + e*/;^(*')'^*'^|i)) . (9) 

The presence of modulates the energy splitting of the Hamiltonian, such that in a frame 
rotating with angular frequency f2, the Bloch vector accumulates a random phase shift relative 
to the y-axis at time t. We characterize /3(t) in the frequency domain by the noise power 
spectrum Splco), defined as the Fourier transform of the two-time correlation function of P{t). 
The influence of the noise term enters the measure of qubit coherence at time r by writing 
W{t) = e^^f'^^ where 

oo 

Xir) = - / ^^Fi.r)dco. (10) 
n J uj^ 



In the expression above, F{lot) is known as a filter-function [221 121] which encapsulates the 
experimental conditions under which qubit coherence is measured. In a Ramsey experiment 
the filter function takes the form 

F(wt) = 4sin2(tjr/2), (11) 

such that we expect an approximate Gaussian decay of coherence as a function of the qubit 
free-precession time, when the small angle approximation is valid. The solid line in Fig. [7^ is an 
envelope function for the decay of qubit coherence based on the construction presented above. 
We plot ^{1 — W{t)) such that complete loss of phase coherence results in a measured value of 
i , commensurate with the decaying contrast of Ramsey fringes. This function incorporates the 
filter function appropriate for a Ramsey experiment and the experimental noise spectrum as 
plotted in Fig. [TJ;. Magnetic field fluctuations in our superconducting magnet were measured 
with a solenoid embedded in the magnet system and a spectrum analyzer. However, the 
vertical scale of this spectrum depends on the inductance of the solenoid which is not precisely 
known, and is hence plotted in arbitrary units. Accordingly, the fitting function contains one 
free parameter, a, which is a multiplicative factor characterizing the strength of the noise. 
With this single free parameter we obtain very good agreement between the form of the 
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envelope function and the measured decay of Ramsey fringes, and find a 1/e decay time for 
fringe contrast T2 w 2.4 ms. 

The experimental noise spectrum displayed in Fig. [Tj: shows that for angular frequencies 
beyond ~2,000 rad/s (f'^SOO Hz), Si3{uj) is reduced by more than four orders of magnitude 
from its low-frequency peak. Therefore, for qubit free-precession times t <C3 ms, the ambient 
magnetic field is approximately constant, and the observed decoherence (recall each point 
represents the results of mutliple experiments that have been averaged together) is due to 
fluctuations in the magnetic field from experiment to experiment (so-called shot-to-shot vari- 
ations). Integrating the measured noise spectrum, and using the scaling parameter extracted 
from fits to Ramsey oscillations, we obtain an rms shot-to-shot variation in the magnetic field 
SB/B ~lxlO-9. 

5.2 Effects of Magnetic Field Inhomogeneities 

Examining the standard deviation of the measurement results for each data point in the 
Ramsey-fringe curves elucidates the dominant mechanisms of decoherence in the system. We 
find that for short times, as described above, the standard deviation oscillates sinusoidally, and 
in a manner consistent with quantum projection noise [45j . As we increase the precession time, 
the sinusoidal variation persists, but the magnitude of the deviation increases commensurate 
with the beginning of decay in fringe contrast. As precession time increases further, we 
observe that the magnitude of the standard deviation begins to decay with what appears to 
be a simple exponential envelope. These measurements suggest that for intermediate times 
(t < 2 ms), qubit coherence is limited by homogeneous dephasing, where all ion spin vectors 
are aligned and the system behaves as an effective single spin. However, as the free-precession 
time is further increased, the effects of sample inhomogeneities become apparent and will 
eventually dominate homogeneous dephasing. 

In our system, inhomogeneity arises from axial and radial gradients in our quantizing 
magnetic field. We characterized the degree of field inhomogeneity by measuring the electron 
spin-flip transition frequency with a small cloud of ions centered at different axial and radial 
positions relative to the fixed trap electrodes. The cloud was moved axially by adding a 
voltage Ay to one endcap and subtracting it from the other. Radial shifts were achieved by 
applying asymmetric voltages to the azimuthally segmented electrodes located in the z — Q 
plane used to create the rotating wall potential. At each location of the cloud a resonance 
profile like that displayed for = tt in Fig. [Jj: was taken, and the electron spin-flip frequency 
was extracted from a fit to the measured Rabi lineshape. 

From these measurements we determined an axial gradient of 9.5 kHz/mm (3.4x10^* 
mT/mm), a maximum gradient transverse to the z-axis of 3.9 kHz/mm (1.4 x 10^* mT/mm), 
and an average quadratic radial dependence of 3.7 kHz/mm^ (1.8 xlO~^ mT/mm'^). For 
single-plane arrays the axial gradient does not produce any inhomogeneous broadening. Oc- 
casionally, however, double-plane arrays with 15 /im spacing were used in experiments, pro- 
ducing a 140 Hz {5B/B ~ 10~^) shift of the qubit resonance frequency between the planes. 
The linear gradient transverse to the magnetic field direction should not produce any static 
frequency shifts due to the averaging effect of the ion array rotation, unlike the quadratic radial 
magnetic field gradient. For a typical array radius of 0.3 mm (corresponding to ~600 ions with 
a 1 kV trap voltage) the quadratic radial field variation produces a 330 Hz {SB/B ^ 2.7x 10^^) 
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shift in the qubit transition frequency between the cloud center and edge. This shift is small 
relative to the linewidth of the qubit transition (for our operational Tjt), but does contribute 
a small rotation error which may dominate our measured operational infidelity. We estimate 
a rotational infidelity of < 1 x 10"'^, consistent with the results of randomized benchmarking 
experiments. 
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Fig. 7. Studies of qubit coherence a) Decay of Ramsey fringes performed with controlled Larmor 
precession. Solid lines are a fit to the envelope of the decay using the coherence integral and filter 
function as introduced in Eg. 1101 and the measured magnetic-field noise presented in panel (c). 
Each data point consists of the average of 20 experiments, b) Standard deviation of averaged 
data in a). Standard deviation oscillates consistent with quantum projection noise. Amplitude of 
standard deviation grows with increasing time, commensurate with a decrease in fringe contrast. 
This indicates the dominance of homogeneous dcphasing processes up to free-precession times of 
~ 1.5 ms. For longer times, magnitude of standard deviation decreases, consistent with the onset 
of inhomogeneous dephasing across the atomic ensemble. Solid lines are guides to the eye showing 
exponential scaling, c) Magnetic field noise measured using a spectrum analyzer connected to a 
shim coil in our NMR magnet, plotted in arbitrary units on a logarithmic scale. Noise power decays 
with an approximate scaling to a low- frequency cutoff near 30 Hz. d) Decay of coherence in 

a spin-locking experiment fit to a simple exponential, and giving an approximate Tip = 688 ms. 
Insets depict pulse sequence and graphical representation of experimental rotation. Rx indicates 
a rotation about the X axis for a variable time. Dashed vertical line indicates the approximate 
experimental duration associated with 200 computational gates in a randomized benchmarking 
experiment, and corresponds to an error ~3x lower than that measured via benchmarking. Z- 
rotations enacted by detuning the microwave frequency for a controlled time while microwave 
switch is off, producing a passive frame change relative to the initial rotational direction. 



5.3 Relaxation in the Rotating Frame 

Homogeneous and inhomogeneous dephasing processes may largely be suppressed by the ap- 
plication of spin echo ,48J and other dynamical decoupling pulse sequences, as will be discussed 
in a later section. However, it is useful to study the fundamental coherence limits of the qubit 
in the absence of extra control pulses, as elucidated by a spin-locking measurement |14[ I49j . 
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The measurement begins with a n/2x pulse which rotates the spin vector to the equatorial 
plane. This is followed by a Tr/2z pulse which rotates the Bloch vector to the X-axis. The 
operation Rx is then performed in which the spin state is continuously driven about the 
X-axis for a variable time S> tv. The experiment concludes with a 3Tr/2z rotation in the 
equatorial plane, and a Tr/2x rotation to the dark state. In this experiment the az operations 
are performed by detuning the microwave drive by 100 kHz for a precise time while the mi- 
crowave switch is off, imparting a controlled Larmor rotation. We observe that error appears 
to accumulate in time as a simple exponential with a decay slow compared to the Ramsey 
free- induction decay time, T2. 

In a spin-locking experiment, the spin vector ideally remains aligned with the applied mi- 
crowave field, but over time the alignment of the Bloch vector with the driving field decays. 
Early Nuclear Magnetic Resonance (NMR) experiments [H] showed that in spin-locking ex- 
periments the net magnetization in the equatorial plane would relax over time. Since the 
spins were aligned with the strong RF driving field, energy relaxation was required to account 
for the change in magnetization, hence leading to the concept of the longitudinal relaxation 
time in the rotating frame, characterized by exponential decay with time constant Tip. In 
these experiments, decay in the rotating frame was attributed to spin-lattice and other similar 
spin-coupling effects. However, due to the large ion spacing in our experiments ('--^lO fim) we 
do not expect such coupling to be significant. Instead, we attribute the observed decay to 
the presence of ambient magnetic field fluctuations. This decay is qualitatively different from 
the effect of magnetic field fluctuations in a Ramsey experiment due to the noncommutative 
nature of the perturbation with the system Hamiltonian in the present setting. While the 
experimental scenario is similar to what is observed in liquid-state NMR, at this time we 
have not developed a complete theory to model the decay of a spin-locking signal in this 
environment. 

5.4 Dynamical Decoupling 

Dynamical decoupling is a quantum control technique designed to extend qubit coherence and 
suppress the accumulation of qubit errors by the repeated application of parity-reversing con- 
trol operations [HI [TH [T71 [HI [TH [201 [211 [25l [26] . A simple and familiar version of dynamical 
decoupling is the Hahn spin echo, widely employed in NMR, electron spin resonance (ESR), 
and quantum information. In a frame co- rotating with the unperturbed qubit splitting, fi, a 
freely precessing qubit in the presence of environmental noise accumulates a random phase 
between the qubit basis states, and hence experiences dephasing relative to the observer. 
However, incorporating a tt pulse halfway through a well-defined free-precession period ap- 
proximately time-reverses the accumulation of phase. So long as fluctuations in the external 
noise are slow relative to the length of the free-precession period, the effect of the noise is 
nearly perfectly cancelled at the end of the free-precession period. This technique was origi- 
nally designed to mitigate the effects of inhomogeneities in NMR experiments |48j . However, 
the description applied above is also appropriate in the case of a single spin. Adding tt pulses 
to the sequence leads to an extension of the coherence time which scales linearly with the 
number of pulses, n 

Considerable research effort has focused on understanding the capabilities of dynamical 
decoupling pulse sequences for suppressing the effects of both dephasing and longitudinal re- 
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laxation, and these techniques are often suggested as candidate error suppression strategies 
allowing for the realization of fault-tolerant quantum computation. Recently, it has been 
discovered that adjusting the relative positions of tt pulses within a multipulse sequence can 
improve the noise suppressing capabilities of a given dynamical decoupling pulse sequence by 
orders of magnitude in particular noise environments [22l [50l [TTl [12] . While dynamical decou- 
pling may be employed to mitigate general decoherence processes, we will focus exclusively 
on dephasing processes in the following discussion. 

The assumption that environmental noise is constant over the duration of the sequence is, 
in general, insufficient, and we must consider how a dynamical decoupling pulse sequence sup- 
presses noise characterized by an arbitrary power spectral density, Sp{oj). We again consider 
a qubit in a noisy environment whose time evolution is given by the Hamiltonian written in 
Eq. [51 As introduced in Eq. [TUl we may write a filter function describing the influence of a dy- 
namical decoupling pulse sequence on qubit coherence [551 [Ml HI] • For any n-pulse sequence 
one may write a time-domain filter function, y„(t), with values ±1, alternating between 
each interpulse free-precession period. The time-domain filter function accounts for phase 
accumulation under the influence of external noise, with the accumulation being effectively 
time-reversed by application of nx pulses (yielding the alternating sign changes). We account 
for nonzero t^^ by inserting delays with this length and value zero into Unit) between free- 
precession periods. Moving to the frequency domain we may then write F{ujt) = |2/„(wt)P, 
where y„(wr) is the Fourier transform of the time-domain filter function, and r is the total 
sequence length. 

Following this construction, for an arbitrary n-pulse sequence with length r we may write 
the filter function in the frequency domain 



F{ujt) = |y„(c^T)|2 = 



1 + (-1)"+^.^'^^ + 2^(-l)^e*'^^'^^ cos (cjT„/2) 



(12) 



where Sjt is the time of the center of the j^^ ttx pulse, and r is the sum of the total free- 
precession time and 7r-pulse times. It is thus apparent that modifying the relative pulse 
locations for arbitrary sequence length leads to modification of the filter function. As F{lut) 
enters the coherence integral, x(''')i ^ multiplicative factor, we see that the spectral overlap 
between Splco) and the filter function determines the coherence of the qubit. Modifying the 
values of Sj therefore provides a mechanism by which the noise suppressing properties of a dy- 
namical decoupling pulse sequence can be tailored to improve qubit coherence. Experimental 
validation of this construction and technique was provided in [11) I12j . 

We present data studying the effects of two distinct pulse sequences as will be described 
in order below: CPMG, after Carr Purcell Meiboom Gill, and UDD, or Uhrig Dynamical 
Decoupling. CPMG is an extension of the Hahn spin echo to a multipulse form, incorporating 
evenly spaced tt pulses about an axis rotated 90 degrees from the direction imparting the initial 
■7t/2x- The UDD sequence is based on Uhrig's discovery [52] that for an n-pulse sequence, 
enforcing certain constraints on the filter function yielded a sequence which improved the 
error suppression characteristics of the sequence in the so-called "high-fidelity" regime (in 
which quantum computations are likely to be performed) by orders of magnitude relative to 
other known sequences. UDD does a particularly good job at suppressing errors associated 
with noise power spectra that have sharp high-frequency cutoffs [22l [23l [5T1 [50l [24l [52l mi [T2] . 



M. J. Biercuk, H. Uys, A. P. VanDevender, N. Shiga, W. M. Itano, and J. J. Bollinger 




Angular Frequency (rad/s) 



Fig. 8. Coherence preservation via dynamical decoupling, (a) CPMG (solid markers) and UDD 
(open markers) sequence performance in ambient noise for various pulse numbers in a Ramsey 
interferometer, as a function of total free precession time. Qubits are rotated to |J.} at the end of 
an experiment in the absence of dephasing. The vertical axis corresponds directly to a normalized 
fluorescence count rate, and hence the probability of completing the experiment in |t) due to 
dephasing. A value of 0.5 corresponds to complete dephasing. Here, = 185 /is. Curves offset 
vertically for clarity. Solid lines represent theoretical fits using the appropriate filter function 
construction as described in Eq. 1121 UDD and CPMG perform similarly for all n in the presence 
of ambient noise with a soft high-frequency cutoff, (b) Similar data taken with an artificially 
engineered noise environment, created using the noise generation system depicted schematically 
in Fig. [3J3. Spectrum is Ohmic (5^(0;) oc uj) with a sharp high-frequency cutoff at 500 Hz. 
In this regime UDD outperforms CPMG, providing significant phase-error suppression over a 
wide range of free-precession times. Error bars on curve for n = 6 derived from normalized 
standard deviation of averaged experiments and demonstrate statistical significance of differences 
in sequence performance, (c) Ambient (dotted line) and Artificial Ohmic (solid line) noise power 
spectra measured directly, and plotted in real units. Ambient spectrum measured using shim coil 
in NMR magnet and normalized using a fit to contrast decay of Ramsey fringes (Fig. [71 a). Ohmic 
spectrum derived from measured frequency noise on PLL reference carrier frequency, acquired 
using a commercial phase- noise detection system. For details see |lll I12| . 
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For these experiments we employ our noise-engineering capabilities and test the error 
suppressing properties of CPMG and UDD in a variety of experimentally relevant noise envi- 
ronments, as demonstrated originally in I12j. In Fig. [8] we show experimental results for 
two different noise environments as well as theoretical fits to our data based on the measure 
of qubit coherence and the filter function introduced in Eqs. [TOl [El These experiments are 
conducted in a manner similar to the Ramsey coherence measurements described in an earlier 
section, where the free-precession period is replaced with a sequence of tt pulses separated 
by free-precession periods. We demonstrate that it is possible to both extend the qubit co- 
herence time by increasing n, but also that the spectral character of Sf}{uj) impacts sequence 
performance as highlighted in Fig. [5^-b. Theoretical fits show good quantitative agreement 
with the data and incorporate a single free parameter — a, the strength of the noise. In 
these fits, a variety of observed nonlinear features are accurately reproduced by theoretical 
expressions incorporating the actual Sf}{u!). For example, in the ambient spectrum, a sharp 
feature near 153 Hz is fully responsible for a rounded plateau appearing in the appropriate 
coherence curves at intermediate times. 

Further studies presented in references pTl [53] also demonstrated our ability to develop 
new suites of locally optimized dynamical decoupling pulse sequences which provide superior 
error suppression relative to all other known pulse sequences. These experiments validate a 
substantial body of literature in the quantum information community and show the versatil- 
ity of our experimental system for the development of novel and powerful quantum control 
techniques. 

6 Nuclear Spin Manipulation 

Experimental demonstrations discussed thus far have focused on high-fidelity quantum control 
using a 124 GHz electron-spin-flip transition. However, this is not the only two-level manifold 
capable of serving as a qubit in this atomic system. ^Be+ has nuclear spin / = 3/2 which 
means that there are three nuclear-spin-flip transitions for each electronic spin state (see 
Fig. [9]) . A state associated with any of these nuclear-spin- flip transitions can be coupled to 
the optically bright state (|t)) for the purposes of detection using multiple rf or microwave 
pulses. Three of these transitions, 

fl : li) = |F = 2,m/ = +3/2, rnj = -1/2) ^ \F ^ 2,m/ = +1/2, mj = -1/2) 
f2 : It) = |F = 2,m/ = +3/2, mj = +1/2) ^ = l,m/ = +1/2, mj = +1/2) (13) 
fS : |F = 2, mi = +1/2, mj = +1/2) ^ = l,m/ = -1/2, mj = +1/2) 

at -288 MHz, -340 MHz, and -286 MHz respectively (at 4.5 T) [54], are either directly 
coupled, or may be coupled using a single pulse, to the optically bright electron-spin-flip 
manifold, allowing for high-fidelity state initialization and projective measurement through 
appropriate mapping between the manifolds. 

The structure of the energy levels, as represented in Fig. [S^ provides these states with 
low-order immunity to magnetic field fluctuations relative to the main electron-spin-flip qubit 
transition, making them attractive for long-term information storage, and similar to those 
used in previous experiments on high-stability frequency clocks |55j . The transitions fl, f2, 
and f3 have been driven using RF fields, and we show a Ramsey experiment on the f2 transition 
(fl and fS not shown) with high contrast for free-precession periods of four seconds. These 
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Fig. 9. Coherent control over a nuclear spin flip in the ground-state manifold of ^Be"*". (a) 
Ground-state energy-level evolution with applied magnetic field. Nuclear spin flip transitions fl, 
f2, and f3 designated in the figure along with the corresponding (m/,mj) quantum numbers, 
fl ~ 340 MHz, f2 ~ 288 MHz, fS ~ 286 MHz. Diagram not to scale, (b) Ramsey fringes measured 
on the f2 nuclear-spin-fiip transition as a function of RF frequency offset from 288,172,932.220 Hz. 
Frequency tuning range of 0.4 Hz is within bandwidth of tt/2 pulses. As one state corresponds to 
the bright state (|t)) in the electron-spin-flip manifold, we directly probe state population through 
the measurement of resonance fluorescence. Total free precession time is 4 s, giving fringe contrast 
'-^90 %. Microwaves are held at a fixed frequency for the duration of the experiment and the RF 
switch provides >60 dB attenuation in the off-state. 
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experiments employ the RF system depicted schematically in Fig.[3t, as well as a loop antenna 
proximal to the vacuum envelope of the trap. 

Unfortunately, the nuclear spin states are weakly coupled to driving fields because of the 
relatively small value of 5/, requiring large RF power to drive rotations at speeds comparable 
to the electron-spin- flip transition. In Ramsey experiments we realize tt/2 times near 0.5 
s, making these operations unreasonably slow relative to the speed of rotations available on 
the electron-spin-flip qubit. However, we believe that careful engineering of our RF antenna 
system will allow the realization of nuclear spin-flip times well below 100 ms. 

7 Future Outlook and Conclusions 

We have presented a set of experiments demonstrating high-fidelity quantum control using 
trapped ions in a Penning trap. Beyond what has been presented in this manuscript and 
in references |1H I12|, I53j . considerable opportunities exist for exploiting the extraordinary 
capabilities of this system in future quantum control and quantum information experiments. 
In this section we provide a brief discussion of future experimental avenues of interest to us. 

7.1 Improved Operational Fidelity 

Effects due to magnetic field inhomogeneity are experimentally observable (see Fig. [7|, al- 
though they could be minimized by working with small ion arrays. The measured magnetic 
field gradients are an order of magnitude larger than those anticipated from the specifications 
issued by the magnet manufacturer. This could be accounted for by the fact that the center 
of our trap is located ~5 mm below the center of the magnet shim coils, and shimming the 
magnet at this point could have produced a much smaller region of field homogeneity. Sig- 
nificantly smaller field gradients might be achieved by moving our trap and re-shimming the 
magnet. This could enable significantly larger arrays to be used in future studies, improving 
operational and measurement fidelity. 

Of primary interest is the ability to augment the capabilities and performance of our 
microwave system. At present our experiments use square microwave pulses implemented 
via TTL control over a p-i-n diode microwave switch. Many interesting quantum control 
experiments may be performed if we have greater control, in real time, over the amplitude of 
the microwave driving field [56l |25l [26| ■ Accordingly, we are interested in the possibility of 
using a voltage-controlled attenuator to allow pulse shaping, and even the possibility of new 
dynamical decoupling techniques that use continuously varied drive fields rather than discrete 
pulses. 

Operational fidelity may be significantly improved by reducing the minimum achievable 
Tjr, currently ~ 90 /is. This is limited by the power output of the Gunn diode and our ability 
to focus the microwaves on the ion crystal. At present our microwave beam waist is ~7 mm at 
the ion crystal, as determined by the teflon lens we employ. We have designed and fabricated 
a new teflon lens that should tighten the microwave beam waist by a factor of two, and hence 
reduce the area by a factor of four. The value of r^r scales linearly with the amplitude of 
the microwave field incident on the ions, and hence a 4x improvement in power density will 
be translated to a twofold reduction in Ttt. It may also be possible to achieve up to a 5x 
increase in microwave power generation by replacing the Gunn diode oscillator used in our 
system. Further, the requirements of our microwave phase-lock-loop limit us to launching 
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only approximately 10 % of the Gunn diode's output from the Horn antenna, a limitation 
that may be lifted by re-engineering our microwave PLL. Overall, we believe it is possible to 
achieve more than a lOx reduction in r^. Finally, if becomes a critical parameter in future 
experiments, reducing the magnetic field will allow access to a wide range of commercial 
microwave components in the W-band near 90 GHz, including higher-ouput Gunn diodes. In 
parallel, timing infidelity issues may be mitigated by employing pulse generators with clock 
frequencies higher than the 20 MHz used through most of our experiments. 

The magnitude of integrated shot-to-shot magnetic field variations represents an infidelity 
of approximately 1 x 10""^ for 200 fj,s. Increasing the length of the magnetic field vector 
due to the driving field in the rotating frame by 10 x will result in a 100 x improvement in 
operational fidelity for fixed magnetic field noise. This will allow us to realize operational 
fidelities well below 1 x lO"'^, appropriate for tests in the "fault-tolerant" regime. 

7.2 Improved Detection Fidelity 

Our present experimental set-up uses an //5 solid angle collection of the ion fluorescence 
scattered perpendicular to the magnetic field. The light collection is limited primarily by the 
size of the hole in the ring electrode which provides optical access. //2 light collection may 
be achieved parallel to the magnetic field, but in the current set-up this access is blocked by 
the teflon lens used for focusing the 124 GHz microwaves on the ion array. A light collection 
of better than f/2 (while simultaneously admitting microwaves) should be possible in a new 
trap, improving our light collection by more than a factor of 6. 

In the absence of dark counts or stray laser light scattered off trap electrodes, detection 
fidelity is typically limited by the mean number of photons scattered by the bright state in 
the amount of time required for the dark state to be repumped to the bright state by the 
detection laser. For the 124 GHz electron spin qubit discussed here, this is the number of 
photons that are scattered on the cooling transition (see Fig.[2|) before the |J,) = (+3/2, —1/2) 
state is repumped to the |t) = (-1-3/2,-1-1/2) state. In general, this scattering rate improves 
with increased magnetic field due to the energy splitting of the excited state in ^Bc+. At 
4.5 T the cooling laser is detuned from the repumping transition by 40 GHz. 

With an //2 light collection and assuming a 20 % quantum efficiency we estimate that on 
average ^20,000 photons can be detected for each bright-state ion before ||) is repumped to 
It). With an array of a few hundred ions this should permit detection infidelities at the 10~^ 
level, allowing the observation of quantum jumps when most of the ions are in the dark state. 

7.3 Exploitation of the Nuclear Spin 

We have directly driven the RF nuclear spin flip transitions discussed in the last section of 
this manuscript, but to date we have not employed them as independent qubit manifolds with 
coherent mapping to and from the electron-spin-flip transition in an algorithmic setting. We 
envision both long-term information shelving [57] , and potentially even multiqubit algorithms 
and advanced quantum control procedures which involve multiqubit interaction [58j . These 
capabilities are predicated on our ability to achieve stronger coupling of the nuclear spins to 
the RF driving field, which will require substantial, but straightforward modifications to the 
RF antenna system. 
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7.4 Single Ion Addressing 

Many quantum information studies can be pursued without addressing individual ions. Exam- 
ples include the dynamical decoupling studies discussed here, and the simulation of quantum 
systems [SHI ED] through the generation of large-scale entanglement. More complex exper- 
iments, including procedures outlined in |61l . may eventually require single-ion addressing. 
While challenging, it may be possible to use strobed or co-rotating laser-beams locked to the 
frequency of the rotating wall potential, similar to the techniques used for single-ion imaging. 
We discuss here the prospects for individual ion addressing with well focussed laser beams 
directed perpendicular to the ion array. 

In Section 13.11 we described the stability of the ion crystal in the frame of the rotating 
wall potential. The periods between "slipping" events — tens of seconds — are sufficiently 
long to make single ion addressing with well focused, rotating laser beams feasible. Crossed 
acousto-optic modulators, combined with an imaging system, can be used to generate a tightly 
focussed laser beam which is displaced from, but rotates about, the rotation axis of the planar 
array. This technique is similar to that demonstrated with stirred Bose-Einstein condensates 

MM- 

Clearly, increasing the length of the time periods where the planar array is locked to the 
rotating wall improves the feasibility of single ion addressing and motivates further studies 
like that discussed in [42]. For instance, the use of a strong rotating quadrupole potential 
which generates a strong perturbation to the ion array boundary should increase the length 
of the "sticking periods" . This regime may be reached with a minimization of the distance 
between the ions and the electrodes imparting the rotating wall potential. It will also be 
useful to reduce the deleterious effects of the perpendicular laser beam torque. This may be 
achieved by appropriate tailoring of the perpendicular laser beam profile; for example use of 
a large beam waist with frequency dispersion across the waist that matches the ion Doppler 
shifts due to the plasma rotation. 

Rather than precisely rotating a well focussed laser beam as discussed above, single ion 
addressing should also be possible with a fixed laser beam and an appropriate center-of-mass 
(COM) shift in the planar array. Consider an ion located a distance tq from the center of a 
single species planar array whose rotation frequency ujr is set by a rotating (e.g. quadrupole) 
potential. The position of the ion can be parameterized by r{t) = j-f^^-^i^rt+s) ^ ^^^^ 
apply a rotating electric field which drives the COM motion of the ions. By rotating the 
field at cOr and setting its strength and phase appropriately, the COM displacement can be 
adjusted to equal —r{t). In this case the designated ion is fixed while the remaining ions 
rotate about it, which allows addressing with a fixed laser beam. This addressing scheme is 
simpler than rapidly rotating a laser beam, but it faces some limitations. First, the COM 
displacement of the array would likely have to be performed slowly relative to the transverse 
mode frequencies of the array, likely limited by the presence of shear modes. These modes 
have not been studied, but their frequencies are reduced by E x B dynamics relative to that 
of an unmagnetized planar array. Also, the presence of impurity ions of heavier mass will 
complicate the COM displacement of the array. A sufficiently large number of impurity ions 
of heavier mass will shield the applied electric field and result in no net displacement of the 
Be+ ions jlO]. 
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7. 5 Conclusion 

In summary, we have demonstrated high-fidelity control over a qubit realized as trapped ^Be+ 
ion crystals in a Penning trap. This trap employs static electric and magnetic fields in order 
to provide charged-particle confinement and yields distinct benefits over RF traps in realizing 
large ion arrays. Qubits are defined using an electron-spin-flip transition at 4.5 T, and are 
driven using a 124 GHz microwave system, providing control over all axes of the Bloch sphere, 
and long coherence times. The average infidelity per randomized computational gate is mea- 
sured through a benchmarking technique to be 8± 1 x 10~^, providing significant advantages 
for studies of high-fidelity quantum control. As an example, we present studies of dynamical 
decoupling pulse sequences, highlighting our ability to suppress the accumulation of phase 
errors and artificially engineer the noise environment in the system. Taken together, these 
experiments illustrate how ion crystals in a Penning trap are a broadly useful model quantum 
system. Future experiments show a path towards higher-fidelity gates, and potential studies 
of large-scale entanglement and quantum simulation. Such experiments are challenging tech- 
nically, but may further improve the value of ion crystals in Penning traps as a fundamental 
tool for studies of quantum information and quantum control. 
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